The Proto-neutron Star Phase of the Collapsar Model and the Route to Long-soft Gamma-ray Bursts and Hypernovae

Recent stellar evolutionary calculations of low-metallicity massive fast-rotating main-sequence stars yield iron cores at collapse endowed with high angular momentum. It is thought that high angular momentum and black hole formation are critical ingredients of the collapsar model of long-soft gamma-ray bursts (GRBs). Here, we present 2D multi-group, flux-limited-diffusion MHD simulations of the collapse, bounce, and immediate post-bounce phases of a 35-Msun collapsar-candidate model of Woosley & Heger. We find that, provided the magneto-rotational instability (MRI) operates in the differentially-rotating surface layers of the millisecond-period neutron star, a magnetically-driven explosion ensues during the proto-neutron star phase, in the form of a baryon-loaded non-relativistic jet, and that a black hole, central to the collapsar model, does not form. Paradoxically, and although much uncertainty surrounds stellar mass loss, angular momentum transport, magnetic fields, and the MRI, current models of chemically homogeneous evolution at low metallicity yield massive stars with iron cores that may have too much angular momentum to avoid a magnetically-driven, hypernova-like, explosion in the immediate post-bounce phase. We surmise that fast rotation in the iron core may inhibit, rather than enable, collapsar formation, which requires a large angular momentum not in the core but above it. Variations in the angular momentum distribution of massive stars at core collapse might explain both the diversity of Type Ic supernovae/hypernovae and their possible association with a GRB. A corollary might be that, rather than the progenitor mass, the angular momentum distribution, through its effect on magnetic field amplification, distinguishes these outcomes.Comment: 5 pages, 1 table, 2 figures, accepted to ApJ

Electromagnetohydrodynamics

Interaction of plasma flow with a magnetic obstacle is a frequent process in many laser-plasma experiments in the laboratory, and is an important event in many astrophysical objects such as X-ray pulsars, AGN, GRB etc. As a result of plasma penetration through the magnetic wall we could expect a formation of magnetohydrodynamic (MHD) shock waves, as well as of electromagnetic (EM) ones. To study these processes we need equations following from hydrodynamic and Maxwell equations, which in the limiting situations describe MHD and EM waves, and are valid for the general case, when both phenomena are present. Here we derive a set of equations following from hydrodynamic and Maxwell equations, without neglecting a displacement current, needed for a formation of EM waves. We find a dispersion equation describing a propagation of a weak linear wave in a magnetized plasma along the $x$ axis, perpendicular to the magnetic field $H_y(x)$, which contains MHD, hydrodynamic and EM waves in the limiting cases, and some new types of behaviour in a general situation. We consider a plasma with zero viscosity and heat conductivity, but with a finite electric conductivity with a scalar coefficient.Comment: 8 papers, 8 figures, 1 table, to be submitted in PR

Cosmic Rays from PeV to ZeV, Stellar Evolution, Supernova Physics and Gamma Ray Bursts

The recent success of a proposal from some time ago to explain the spectrum of cosmic rays allows some strong conclusions to be made on the physics of supernovae: In the context of this specific proposal to explain the origin of cosmic rays, the mechanism for exploding supernovae of high mass has to be the one proposed by Bisnovatyi-Kogan more than 30 years ago, which was then based on a broader suggestion by Kardashev: A combination of the effects of rotation and magnetic fields explodes the star. Interestingly, this step then leads inevitably to some further suggestions, useful perhaps for the study of gamma ray bursts and the search of a bright standard candle in cosmology

A Two-Dimensional MagnetoHydrodynamics Scheme for General Unstructured Grids

We report a new finite-difference scheme for two-dimensional magnetohydrodynamics (MHD) simulations, with and without rotation, in unstructured grids with quadrilateral cells. The new scheme is implemented within the code VULCAN/2D, which already includes radiation-hydrodynamics in various approximations and can be used with arbitrarily moving meshes (ALE). The MHD scheme, which consists of cell-centered magnetic field variables, preserves the nodal finite difference representation of div(\bB) by construction, and therefore any initially divergence-free field remains divergence-free through the simulation. In this paper, we describe the new scheme in detail and present comparisons of VULCAN/2D results with those of the code ZEUS/2D for several one-dimensional and two-dimensional test problems. The code now enables two-dimensional simulations of the collapse and explosion of the rotating, magnetic cores of massive stars. Moreover, it can be used to simulate the very wide variety of astrophysical problems for which multi-D radiation-magnetohydrodynamics (RMHD) is relevant.Comment: 22 pages, including 11 figures; Accepted to the Astrophysical Journal. Higher resolution figures available at http://zenith.as.arizona.edu/~burrows/mhd-code

Features of the Acoustic Mechanism of Core-Collapse Supernova Explosions

In the context of 2D, axisymmetric, multi-group, radiation/hydrodynamic simulations of core-collapse supernovae over the full 180$^{\circ}$ domain, we present an exploration of the progenitor dependence of the acoustic mechanism of explosion. All progenitor models we have tested with our Newtonian code explode. We investigate the roles of the Standing-Accretion-Shock-Instability (SASI), the excitation of core g-modes, the generation of core acoustic power, the ejection of matter with r-process potential, the wind-like character of the explosion, and the fundamental anisotropy of the blasts. We find that the breaking of spherical symmetry is central to the supernova phenomenon and the blasts, when top-bottom asymmetric, are self-collimating. We see indications that the initial explosion energies are larger for the more massive progenitors, and smaller for the less massive progenitors, and that the neutrino contribution to the explosion energy may be an increasing function of progenitor mass. The degree of explosion asymmetry we obtain is completely consistent with that inferred from the polarization measurements of Type Ic supernovae. Furthermore, we calculate for the first time the magnitude and sign of the net impulse on the core due to anisotropic neutrino emission and suggest that hydrodynamic and neutrino recoils in the context of our asymmetric explosions afford a natural mechanism for observed pulsar proper motions. [abridged]Comment: Accepted to the Astrophysical Journal, 23 pages in emulateapj format, including 12 figure

A supersymmetric model for triggering Supernova Ia in isolated white dwarfs

We propose a model for supernovae Ia explosions based on a phase transition to a supersymmetric state which becomes the active trigger for the deflagration starting the explosion in an isolated sub-Chandrasekhar white dwarf star. With two free parameters we fit the rate and several properties of type Ia supernovae and address the gap in the supermassive black hole mass distribution. One parameter is a critical density fit to about $3 \cdot 10^7$ g/cc while the other has the units of a space time volume and is found to be of order $0.05\,$ Gyr $R_E^3$ where $R_E$ is the earth radius. The model involves a phase transition to an exact supersymmetry in a small core of a dense star.Comment: 20 pages, 5 figures, expanded version to be published in Physical Review

Application of ESI FT-ICR MS to Study Kraft Lignin Modification by the Exoenzymes of the White Rot Basidiomycete Fungus TrametesHirsutaLE-BIN 072

Trameteshirsuta is a wood rotting fungus that possesses a vast array of lignin degrading enzymes, including7 laccases, 7 ligninolyticmanganese peroxidases, 9 lignin peroxidases and 2 versatile peroxidases. In this study,electrospray ionization Fourier transform ion cyclotron resonance mass spectrometry (ESI FT-ICR MS)was used to examine kraft lignin modification by the enzymatic system of this fungus.The observed pattern of lignin modification suggested that before the 6th day of cultivation,the fungal enzymatic system tended to degrade more oxidized moleculesand, hence, less recalcitrant molecules, with the production of hard-to-modify reduced molecular species. At some point after the 6th day of cultivation,the fungal enzymatic system tended to degrade more oxidized moleculesand, hence, less recalcitrant molecules, with the production of hard-to-modify reduced molecular species. At some point after the 6th day of cultivation,the fungus started to degrade less oxidized, more recalcitrant, compounds, converting them into the more oxidized forms. The altered pattern of lignin modification enabled changes in the fungal enzymatic system. These changes were further attributed to the appearance of the particular ligninolyticmanganese peroxides enzyme(MnP7), which was added by the fungus to the mixture of enzymes that had already been secreted (VP2 and MnP5). Keywords: wood rotting fungi, kraft lignin, mass spectrometry, peroxidase

Equilibrium Configurations of Strongly Magnetized Neutron Stars with Realistic Equations of State

We investigate equilibrium sequences of magnetized rotating stars with four kinds of realistic equations of state (EOSs) of SLy (Douchin et al.), FPS (Pandharipande et al.), Shen (Shen et al.), and LS (Lattimer & Swesty). Employing the Tomimura-Eriguchi scheme to construct the equilibrium configurations. we study the basic physical properties of the sequences in the framework of Newton gravity. In addition we newly take into account a general relativistic effect to the magnetized rotating configurations. With these computations, we find that the properties of the Newtonian magnetized stars, e.g., structure of magnetic field, highly depends on the EOSs. The toroidal magnetic fields concentrate rather near the surface for Shen and LS EOSs than those for SLy and FPS EOSs. The poloidal fields are also affected by the toroidal configurations. Paying attention to the stiffness of the EOSs, we analyze this tendency in detail. In the general relativistic stars, we find that the difference due to the EOSs becomes small because all the employed EOSs become sufficiently stiff for the large maximum density, typically greater than $10^{15}\rm{g} \rm{cm}^{-3}$. The maximum baryon mass of the magnetized stars with axis ratio $q\sim 0.7$ increases about up to twenty percents for that of spherical stars. We furthermore compute equilibrium sequences at finite temperature, which should serve as an initial condition for the hydrodynamic study of newly-born magnetars. Our results suggest that we may obtain information about the EOSs from the observation of the masses of magnetars.Comment: submitted to MNRA